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ABSTRACT 


The hurricane surge problen and the relative importance of the 
parameters involved are discussed. Three practical methods of forecast- 
ing the maximum height of the surge are explained, and tested for reli- 
ability on nine hurricanes covering the period 1954 through 1959. A 
tentative sclution is presented for making a fourth methed of height 
forecasting operational, and its merits relative to the first three 
methods are investigated. A quantitative method of forecasting the 
time of occurrence of the maximum surge is developed. This method is 
testec and the results are discussed. 

The writer wishes to express his appreciation for the assistance 
and encouragement given him by Professor Glenn H. Jung of the U. S. 
Naval Postgraduate School in this study. 

The author is also indebted to Mr. D. Lee Harris of the U. S. 
Weather Bureau for his suggestions and invaluable aid in supplying 


needed data upon which parts of this study were based. 
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1, BACKGROUND 


ieee eniuroduc t10nN . 

The rise in sea level associated with the passage of a hurricane 
accounts for the largest portion of all the damage and loss of life 
attributed to such storms nee Cr 11] « The relative importance of 
this part of the overall hurricane problem has been known for some time, 
but it has only been within the past five or six years, since the very 
destructive Atlantic coast storms of 1954 and 1955, that intensive 
studies have been conducted in the United States to obtain practical 
methods of predicting the behavior of the sea along the coast during 
mie passage ol a hurricane. 
ee Local Bilecus,. 

The problem of deriving any reasonably accurate prediction methods 
is greatly complicated by the very important effects of local topography, 
both above and below the mean water level. These local effects take the 
form of extra high water at the heads of converging bays and at exposed 
points, lower levels where bays and estuaries are wide oa » and some=- 
times violent oscillations of the water levels along the coast for long 
periods after the passage of a storm. 

If these local effects”are considered in™a forecasting technique, 
then the area over which the method can be used is greatly reduced. 

Some "small-area"! techniques have been developed, such as Wilson's 
technique for New York Bay fala) 3 however, due to the limited applica~ 
tion of such methods, area-wise, they were nct considered in this study. 

Most prediction techniques have attempted to eliminate local effects 
as much as possible and can therefore be used over long stretches of 


coastline. For the purposes of this paper only these "wide-area" tech= 
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niques will be discussed. 
1.3 Coastal Water Levels. 

There are three fairly distinct reactions of the sea level as the 
result of a hurricane which crosses the coast. These are shown in 
Figure 1. 

The forerunner is a very gradual rise in sea level which occurs a 
few days before the arrival of a hurricane, when it is still more than 
200 miles from the coast. This build-up can, under certain conditions, 
cause flooding of a low-lying coast; however, this is not normal. 

Following the forerunner, and within several hours of the passage 
of the hurricane, there is a very rapid rise in the water level. This 
is the storm surge which can be as high as eight to ten feet above the 
normal water level. 

After the storm surge reaches its peak the water drops rapidly 
towards its normal level, but often "over-shoots" this mark; the water 
level hits a low point, builds up again, and continues to oscillate 
for a period of time. This is the resurgence which appears to be due 
to the free motion of water as it returns to its normal level [10 | ° 
In closed bays and estuaries this oscillation of the water can become 
even higher than the storm surge, and under such local conditions it 
can become a very serious problem; however, the resurgence is usually 
of small magnitude along the open coast. 

Since the forerunner and the resurgence normally are not very 
destructive, the present state of prediction of high water levels is 
generally confined to the storm surge. Accordingly the remainder of 
this paper will deal only with the desired parameters relating to the 


storm surge. 
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1.4 Storm Surge Definition, 


The meaning of "storm surge" used in connection with this study 
is given by the definition of the National Hurricane Research Project 
of the U.S. Weather Bureau 16 | » in this interpretation the storm 
surge is the actual water level less the astronomical tide and secular 
anomolies in the sea level. Figure 2 shows the effects of these three 
components on the actual water level for a general case. 

The astronomical tide is the familiar diurnal or semi-diurnal 
change in sea level which can be obtained from tide tables. Generally 
the range of the astronomical tide is greater on the Atlantic coast than 
on the Gulf coast; however, it is of great importance in both cases due 
to the lower land elevations in the Gulf area. The astronomical tide 
becomes an important factor through its phase relationship with the 
storm surge. If the surge occurs when the astronomical tide is low, 
the actual water level may not be high enough to cause much harm; how- 
ever, if the surge occurs at high tide, extremely high water levels may 
result with correspondingly large damage along the coast. If the fore= 
caster knows when to expect the storm surge component, he can determine 
its phase relationship with the astronomical tide. 

In the absence of a storm the secular anomoly is the difference 
between the astronomical tide and the actual tide, or essentially the 
level of the water that the tide table does not account for. It can 
amount to a foot or more in height and usually lasts for several weeks 
at a time. The present method for determining this anomoly requires 
the forecaster to be familiar with the predicted and actual tides for 
a period of a week preceeding the forecast time. He then makes a 


qualitative estimate of the anomoly from the differences between these 
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two quantities. A program has recently been established by the U. S. 
Weather Bureau in cooperation with the Coast and Geodetic Survey where= 
by remote recording tide gages have been installed at many Atlantic and 
Gulf Coast weather stations. These stations are provided with predic= 
tions of the astronomical tide and in turn they transmit the difference 
between the actual and predicted values along with their synoptic reports. 
When the program is fully operational a central office will apply an 
averaging technique to determine the final secular anomoly, which will 
then be available for more closely predicting actual tide levels, 

The final, and most important component, is the storm surge itself. 
This is the change in the sea level that is the direct result of the 
hurricane high winds, low pressures, and the additional piling up of 
water near the shore by breaking waves. 

1,5 Storm Surge Parameters. 

In order to promulgate a satisfactory warning with regard to the 
high waters associated with a hurricane, it would be necessary for the 
forecaster to describe the distribution of water levels along the coast, 
the time when the water would reach its maximum height, and the duration 
of the high water. This problem can be broken down into the forecasting 
of five different parameters involving the storm surge. 

These are: 

lL. The maximum height. 

2. The location of the maximum height. 

3. The lateral distribution of the surge along the coast. 

4. The time that the maximum surge occurs. 

5. The duration of the surge. 


Figure 3 shows these parameters with values for a very general "average" 
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2ase. 

The most important surge parameter is the maximum height that it 
will reach, gince the other parameters will have little meaning other-= 
wise. This is the only parameter quantitatively predictable at the 
present time, 

Once that the maximum surge has been predicted, the forecaster 
would next be interested in knowing the time that the maximum would 
occur. A prediction of this parameter would not only allow timely 
action to be taken to protect lives and property; but, even more 
important, it would permit correct phasing of the surge with the 
astronomical tide so that the actual water level could be obtained 
more accurately. At present there are only qualitative rules to guide 
the forecaster in the prediction of the time of the maximum surge. 

The forecaster would also be interested in the remaining three 
parameters, and again he would have to turn to qualitative rules. This 
study has not considered these three parameters. 

1.6 Purpose. 

The remainder of this paper will deal with the two parameters, the 
neight and the time of the maximum surge. The results of a comparative 
study of four different height forecasting methods will be given in 
section 2, and Section 3 will deal with the formulation of a technique 


for quantitatively forecasting the time of maximum surge, 
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< EVALUATION OF HEIGHT FORECASTING METHODS 


2,1 Methods of Approac 


De 





The problem of forecasting the height of the maximum surge result-= 
ing from the passage of a cyclonic storm has been studied for some time 


throughout the world. For many years only a strictly dynamic approach 


> 


conld be made since actual hurricane data were scarce and mostly un= 


¢ 


Other than studies made in the United States, notable contributions 
vc the dynamics of the storm surge problem have been made by Dutch, 
English, and Japanese writers. At present the application of strictly 
dynamical theory to the prediction of the surge associated with storms 
alfecting the United States can not be made with reliability. However, 
the dynamics of the storm surge are still being studied and in the 
ture this theoretical approach may be of great value to practical 
prediction of surge height. 

Jr. recent years there has been considerable interest in collection 
and analysis of hurricane data which could be of use in the development 
of empirical surge forecasting techniques, With these data scme 
progress has been made toward the development of practical methods of 
forecasting the maximum height of the surge. The present methods avail- 
sole are all based on theoretical principles, but are empirically derived. 
z2o2 Methods Considered. 

Four prediction techniques were considered to be evaluated in this 
study. The derivation and use of each of these methods is presented in 
Appendices I through IV. 

The basic ccnsideration in each of these methods is that the height 


of the maximum surge is proportional to the pressure gradient. in the storm, 
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The first two methods use a pressure parameter alone, and the other 
methods supplement this with additional parameters. 

The first practical forecast method to be published was developed 
by Conner, Kraft, and Harris in 1957 12 | - This method will be referred 
to as the "CKH method”. Just a few months after the CKH method had been 
presented, Hoover published a method which was quite similar to the first, 
but based on additional data which had been analysed in a different manner 
| In August of 1959 the U.S. Weather Bureau officially published 
"An Interim Hurricane Storm Surge Forecasting Guide" as one of the re- 
ports of the National Hurricane Research Project [6] « This forecasting 
technique, which will be referred to as the "NHRP method", adds a para- 
meter which takes into account the depth of water over which the hurri- 
cane travels. The final technique that was considered in this study 
was published in December of 1959. It was developed by Kajiura of the 
Agricultural and Mechanical College of Texas EI e in this method two 
parameters are added to the basic pressure parameter. One of these is 
a constant that depends on the location of the station where the maximum 
surge occurs. The second additional parameter introduced in this method 
is the dynamic amplification factor. Since the method, as presented, 
does not give any rules for picking a value of the dynamic amplification 
factor in a particular storm situation, the first investigation in this 
study tested each storm for three values over the possible range of 
this factor as indicated by Kajiura. 
2.3 Data Used. 

In studying the various prediction techniques it was noted that 
different authors had assigned different numerical values to the same 


parameter for the same storm. This is probably due to the fact that as 
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hurricane data have been continually processed and analysed, the values 
which are considered correct naturally have changed from time to time. 
The data used in this study for testing the various methods of surge 
height prediction were obtained from the Climatological Data Periodicals 

[22 | and a series of Storm Surge Data Charts provided by Mr. D. Lee 
Harris of the U. S. Weather Bureau, Some data for Hurricane AUDREY 
were obtained from the National Hurricane Research Project Report 
foe 23 [4]. 

Nine hurricanes covering the period from 1954 through 1959 were 

considered in this study for testing the prediction techniques. These 
storms, the maximum height of the surge, and the location of the maximum 


surge are listed in Table I. 


2.4 Comparison of CKH, Hoover, and NHRP Methods. 


The Kajiura method at present can not be considered as a practical 


forecasting technique since it still has not been determined how the 
on=the-spot forecaster would determine the dynamic amplification factor 
which is required in this method. For this reason the first three 


methods will be compared separately at first, while the Kajiura method 


results enter the comparison later. 


The errors which resulted from applying the CKH, Hoover, and NHRP 


forecasting methods to the nine test storms are shown in Table I. The 
error is taken as the forecast height minus the actual height, so that 


a negative error indicates that the forecast was too small. Table II 
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shows these results tabulated according to how well each method forecast 
for each storm. 

The NHRP method clearly seems to be the most reliable of these 
three forecasting techniques. Although both the CKH and the Hoover 
methods forecast the same number (3) of storms with minimum error, the 
CKH technique seems to be the more consistently reliable, achieving 
second place in the forecast of five storms. 

Although, for certain storms, noteably HAZEL, the error is quite 
large, it should be noted (Table I) that all three of these methods fore~ 
cast nearly the same height of maximum surge for each storm. The larg~ 
est difference for any one storm was only a foot and a half, for hurri- 
cane AUDREY. 

2.5 Kajiura Method, 

As a first test, the Kajiura method was applied to each storm using 
three values of the dynamic amplification factor, r. The errors which re- 
sulted from this test are shown in Table III. When compared with the re= 
sults obtained from the first three methods (Table I), this prediction 
technique appears to be much better. 

An investigation of the dynamic amplification factor was conducted 
in order to determine whether the good results obtained from the Kajiura 
method were due solely to the wide range of heights covered by using 
values of r from 1,0 to 2.0. The dynamic amplification factor has been 
determined for a theoretical one=dimensional storm model by Reid [ 10 | 
and graphed by Kajiura. This graph is reproduced in Figure 4. Kajiura 
has pointed out that the values of r for a two-dimensional model will 
be a little smaller than those shown in the graph. 


In order to enter the graph of the dynamic amplification factor, 
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lst PLACE 3rd PLACE 
FORECAST (fest. with 2nd PLACE (fost ® wath 
METHOD min, error) max, error) 


CAROL EDNA AUDREY 
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ome 
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the width of the shelf and the scale of the storm have to be known, 

In this study the width of the shelf was taken as the distance that the 
hurricane travelled over the continental shelf. The scale of the storm 
1s essentially the fetch of the hurricane winds. As an approximation 

to the fetch the average radius of the hurricane out to the 1000 mb 
isoline was chosen. Using these values, as shown in Table IV, the ampli- 
fication factors for the one-dimensional casewere determined from the 
raph of Figure 4. 

Next a hindcasting technique was applied to the Kajiura formula for 
forecasting the height to determine the actual value of r for each case, 
The actual and one-dimensional values of the amplification factor were 
next compared and it was determined that the one-dimensional values were 
O.8 units too large in the mode. Using the same curve as presented for 
the one-dimensional case, the r scale was shifted 0.8 units as a rough 
approximation to the dynamic amplification factor for the actual storm 
CASE. 

Using the semie=cnsT 0s graph, r was determined, and a forecast of 
the maximum surge height was made for each of the seven applicable storms. 
The errors that resulted from forecasts made using the scale=-shifted graph 
are tabulated along with the errors that resulted from the first three 
methods, in Table V. These results obtained from using the Kajiura 
method with the scale=shifted graph for amplification factor are consider- 
ed along with the results of the other three methods in Table VI, which 
shows how well each of the four methods predicted the height for each of 
the seven Atlantic coast hurricanes. These tables would seem to indicate 


that the Kajiura method is the best of those tested. 
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2.6 Height Forecasting Conclusions. 


It is obvious that the final test of the Kajiura method is far from 
peing unbiased since the selection of the dynamic amplification factor, 
r, was based on the actual maximum surge height for the storms considered. 
The results can not leac to a definite conclusion that the Kajiura method 
is more reliable than the cthers. However, it is believed that further 
study of this method would be well worth-while. The scale-shifted graph 
for determination of r that was used in this study is at best only a crude 
approximation to the actual storm case, since the only change from the one- 
dimensional model curve is the shift of the r scale, and this is based on 
data from only seven hurricanes. An empirical study of this factor us- 
ing more data probably would lead to an entirely different curve, It 
should pe ened out also that the measurement of the "scale of the 
storm" is extremely difficult. Further study should be made in this 
area with an aim toward cefining this measurement or toward substituting 
a different parameter which would give desirable results. 

With additional study, the Kajiura method could become an operational 
forecasting method which may be more accurate than the existing practical 
techniques. However, at the present time, it is believed the NHRP method 


is the most reliable for forecasting the height of the maximum surge. 


eZ 





3. <A QUANTITATIVE TIME FORECASTING METHOD 


3.1 Importance of Time Forecasting. 


Once the practical forecaster has predicted the maximum height of 
the surge, he then would be interested in the time that the maximum 
surge would occur. As pointed out before, a knowledge of the time of 
the maximum surge would aliow more opportune measures to be put into 
effect for the protection of life and property; it also would help the 
forecaster to correctly predict the phase of the surge with relation 
mo Lhe astronomical tide, therefore giving more nearly correct predic- 
tions of the actual water level. This study, therefore, has attempted 
to develop a simple, quantitative technique for forecasting the time 
of the maximum surge. 
pee Development of the Technique. 

For determining the time that the maximum surge would occur, it was 
first assumed that the surge was the result only of the wind-generated 
waves produced at the storm, and that the surge travelled at the group 
velocity of these waves. The basic procedure was to determine an effec- 
tive fetch at the position of the hurricane; permit decay of these waves 
enroute to the coast and determine the resulting period after decay; and 
from this period, calculate the time that would elapse during travel 
over the decay distance. 

The general technique suggested by Francis [3] was followed in 
developing a convenient measurement of the fetch. A simple circular 
hurricane model was adopted as shown in Figure 5. This model was divided 
into eight sectors, each of which can be considered as a fetch which 
emits waves in a direction depending on the orientaticn of the sector 


: ; O 
and the direction of movement of the storm, The 90 sector was chosen 
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as the one from which the highest waves would be emitted, and the one 
that would have the most effect on the storm surge. This is in agree- 
ment with actual wave reports from numerous storms, and with the fact 
that the storm surge is normally to the right of the hurricane path. 
The fetch, at distance Rt from the center of the storm is given by 
Francis [3 | as: 


ee = sin Mp 


wheres A = angular width of sectors 
(A = 45> for model chosen) 


Francis has found, using Sverdrup-Munk relationships, that for wind 
speeds in excess of 40 knots the significant period of the waves is 
essentially the same for both the case of stationary and moving storms. 
(i.e. the same period is produced for all winds above 40 knots blowing 
over the same-sized fetch.) Therefore it was assumed in the develop- 
ment of this technique that the wind speed in the fetch was 40 knots. 
In the actual forecasting case it would not be possible to tell where 
the wind speeds dropped below 40 knots, since seldom are there more 
than three or four simultaneous reports within the hurricane area. 
Based on some detailed wind field analysis conducted by Redfield and 
Miller lo] the area between 1000 and 1005 mb was chosen as the cut- 
off point for 40 knot winds, Both of these values were tested, using 
Wilson's C14 | Ht-FT diagram (Figure 6) for several actual storms 
and the resulting periods were nearly identical. Since the 1000 mb 
isoline almost always appears in map analyses, this was selected as 


the basis for determining R, the radius of the storm in the calicula- 


tPrancis uses the symbol "r", "R'" has been substituted in this 
study for the radius of the storm to avoid confusion with Kajiura's 
ies the dynamic amplification factor. 
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Zion of the fetch, Fee 
The model fetch which will give the necessary information with 
regard to the surge lies 90° to the right of the storm path. Over 
this fetch the wind speed is 40 knots, and the size of the fetch is 
determined by the distance from the center of the storm to the 1000 mb 
isoline, 
The next problem was to show the relationships between the radius 


of the storm out to the 1000 mb isoline, R the speed of movement of 


OO? 
the storn, Cos and the resulting significant periods, Tye 

A series of effective fetches was chosen to cover the possible 
range in the hurricane situation. The corresponding Tp for a 40 knot 
wind blowing over the chosen effective fetches were determined from 
Wilson's Ht-FT diagram | 14 | which is inserted in this report as 
Figure 6, 

The effective fetch, F', equals the true length of the fetch, F, 
plus the distance that the fetch has travelled due to the movement of 
the storm. The distance that the fetch has travelled is equal to the 


minimum duration, t times the speed of the storm. The resulting 


‘min? 
equation for the effective fetch is [3] : 
\ 

(2) pai a C Te 
This equation can be rearranged for the determination of F, knowing 
Ft, t'_.,s and C,. The minimum duration was determined from the Ht-FT 
diagram; and several values of the speed of the storn, Cs were chosen, 
covering the range from 0 to 35 knots. With these values and the pre- 
viously-chosen values of F' the fetches were calculated. The fetch, F, 
was taken as F, in equation (1) and the associated values of R were cal- 


culated. A graph of R vs. Tp for the various storm speeds was prepared. 
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This graph showed that there was not a large range of periods associated 
with the range in storm speeds. Since it was desirable to eliminate as 
many parameters as possible in order to simplify the forecasting proce- 
00 VS« Tp was 


adopted to cover all speeds of storm movement. This curve is presented 


dure, a mean speed was chosen and the associated curve of R 


in Figure 7. These significant periods then had to be modified to account 
for decay of the waves from the end of the fetch into the coastline. 

The Bretschnieder {2 | wave decay curves (Figure 8) were applied to 
varicus values of Te for decay distances from 25 to 1000 n.mi. and the 


decay periods, Ty were calculated. The resulting values were drawn as 


a graph of R,. vs. T,, for various decay distances, D. This graph was 


00 D 


then transformed into nomogram form, 

A second nomogram was drawn for the travel time of waves, which 
depends on the group velocity associated with the period at the end of 
the decay distance [1] » from the formula: 


~ cf Vv 
ey Pelee /y 


where: ty = travel time 


D = decay distance 

Tp = period at end of decay 
Since these two nomograms had two common parameters, Ty and D, they 
could be combined easily. The combined t-D-Th-Roo nomogram is shown 
as Figure 9, If the assumption can be made that over the entire decay 
distance the waves travel in deep water, then the forecaster would not 


be interested in the value of the period, T This assumption is parti- 


i” 
cularly valid when the travel time in shallow water is small compared 
with the total. Therefore, another nomogram was drawn which is the 

same as that shown in Figure 9 except that TH has been eliminated and 


the scalé changed somewhat. This nomogram shown in Figure 10 is the 
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final suggested forecasting diagram. 
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2 "Deep Water’ Time Forecasting Procedure, 





To make a forecast of the time that the maximum hurricane surge 
will occur, the forecaster would proceed as follows: 

Using the latest synoptic hurricane plot; 

1. Measure the distance, in degrees of latitude, 
from the center of the hurricane to the 1000-mb 
isoline in a direction 90° to the right of the fore- 
cast storm track. This is Roos 

2, Measure the distance, in degrees of latitude, from 
the position of the storm in the general direction 
of the forecast storm track to the coast line. 
Thias is D. 

3, Enter the Ryo7tg-D nomogram (Figure 10) with the 
two values measured in steps l. and 2. to obtain 
the time, tae 

4.o Add ty to the time of the synoptic plot of the hurri- 
cane, to obtain the forecast time of the maximum surge. 

304 "Deep-Shallow Water" Time Forecasting Procedure. 

For the purposes of this study the usefulness of making a correc- 
tion for wave travel close to the coast was investigated. To do this 
the period, Tp? was determined from the by~D-TH Roo nomogram (Figure 9) 
and the depth at which the waves would "feel bottom" was calculated from 
the formulas: 


ea (Ty 


© 


where: L deep water wave length (ft) 


1+ % i. 


Q 
i! 


where; depth that waves feel 


bottom 








Figure 10, 
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This depth was located on a chart along the storm path. From the posi- 
tion of the hurricane to this point the waves were considered to trevel 
over deep water and the time was obtained directly from the nomogran, 
For the remaining distance, from the point where the waves feel bottom 
to the coastline, the waves were considered to travel in shallow water 
(i.e., as long waves). The time for this distance was determined by 


doubling the corresponding "deep-water" time, t,, obtained from the nomo- 


d 
gram. 

The time for deep water travel and the time for shallow water travel 
were added and then the forecast time was obtained as previously describ- 
ed. 

A correction for wave travel in intermediate water was not consider- 
ed since one of the objectives in the development of this forecasting 
procedure was to keep it as simple as possible. It is seen that simply 
introducing a correction which assumes a sudden change from deep to 
shallow water has complicated the procedure and requires more "equipment"! 
than before. 

3.5 Evaluation of Time Forecasting Methods. 

The Weather Bureau's Interim Forecasting Guide 6 | Paves amoer.es 
of qualitative rules for prediction of the time of the maximum surge, 
based on the path of the hurricane relative to the coastline. All of 
the storms considered in this study fall into a catagory for which these 
rules indicate that the maximum surge occurs "within a few hours" of the 
closest passage of the storm. In this study, "within a few hours!" was 
interpreted to be plus or minus two hours. The time of the closest 
passage of the storm was taken as the time of the lowest pressure as 


listed in the Climatological Data Periodicals [12 | for the stations 
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where the maximum surge occurred. The results obtained from applying 
this qualitative technique are listed in Table VII. 

Forecasts were also made considering that the surge travelled at 
the same speed as the storm. In this case the surge would arrive at 
the same time as the hurricane. Predictions of the time of the maximum 
surge based on this assumption are also listed in Table VII. The times 
that the hurricanes crossed the coast were obtained from the Storm Surge 
Data Charts, provided by Mr. D. Lee Harris, the Climatological Data 
National Summary [12 | for Hurricane GRACIE, and the National Hurricane 
Research Project Report No. 23 [| for Hurricane AUDREY. 

Both the "deep water" and the "deep=-shallow water" methods, pre- 
viously described, were used to forecast the time of the maximum surge 
for the same nine hurricanes as previously used in the test of the 
height forecasting methods. 

The storm sizes, R,9, and the storm tracks for 1954, 1955, and 
1956 hurricanes were obtained from the Annual Tropical Storm Reports 
of the U. S. Navy Fleet Weather Central [13]. For 1957 and 1959 
hurricanes the storm sizes were measured from the Facsimile Charts 
which had been transmitted by the U. S. Weather Bureau. Storm tracks 
for these latter hurricanes were obtained from the Climatological Data 
Periodicals [12 | - the time of the maximum surge, for all storms 
except AUDREY and GRACIE, were obtained from the Storm Surge Data Charts 
provided by Mr. Harris. This last parameter for hurricane AUDREY was 
obtained from the National Hurricane Research Project Report No. 23 

Le The time of the maximum surge for hurricane GRACIE was taken from 
the September 1959 Climatological Data National Summary tinie) : 


For each hurricane considered,a time of the maximum surge was fore= 
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a a eee ee 

CAROL 31-1040 31-0840 to 1210 
11-1353 11-1153 to 555 
15-1150 15-9959 to 1350 15-1200 
DIANE 17-1600 17-1030, _ 17-0830 to 30 AV~1030 | 
IONE 19-1000 1S-0615 19-0415 to 0215 19-0600 
GRACIE 29-1100 29-1316 29-1116 tO BSG 29-1330 
MeOSSY 24-0330 2h-O1LL5 2372315 to 24-0345 24-0200 
‘aubrey | 27-1100 27-0830 27-0630 tO 1930 2(-0630 


a 





a at 


t - actuel time of the maximum surge 


bo -- time of passage of lowest pressure at station 
of maximum surge 

ual - "qualitative" forecast time of maximum surge 

ve - "Surge speed = storm speed" forecast time of 


maximum surge 


Tadle VII. NQUALITATIVS" & "SURGE SPEED = STCRII SPSED" 
FORECASTS OF TIME OF MAXIMUM SURGE 
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cast for every reported position within 1200 n.mi. of the coast line, 
for which data were available. The results of these forecasts using 
the "deep water" and "deep-shallow water" methods are tabulated in 
Tables VIII & IX. It can be seen that for most of the storms tested 
using the "deep water" method the forecast time approached the actual 
time from the early side with very few "late" forecasts. Using the 
"deep=-shallow water" technique the approach of the forecasted times 
was the same, but many more "late" forecasts resulted when the storm 
was close to the coast. 

The results obtained from all four of these techniques for forecast- 
ing the time of the maximum surge ("qualitative", "surge speed = storm 
speed", "deep water", and "deep-shallow water" methods) can be examined 
by looking at the errors that resulted from the application of each of 
them to the nine storms which were considered. These errors are snovm 
in Table X. The error is taken here as the forecast time minus the 
actual time of the maximum surges; therefore a plus error indicates a 
late forecast, 

Forecasts made using the qualitative rules bracketed the actual 
time in five out of the nine cases. Of course each of these forecasts 
cover a range of four hours which is somewhat undesirable. A compari- 
son of the three remaining techniques, which give a single time for 
each forecast, show that the closest forecasts for more of the hurricanes 
(four out of nine) are obtained from assuming that the surge travels with 
the hurricane. The "deep water" method made the best forecasts for three 
of the storms, and the "deep-shallow water" method was closest for only 
two hurricanes. It should be pointed out that the forecasts made using 


the first two methods were based on actual hurricane data. If forecast 


D4 





time of R D it t (fest- 
max. surge actual) | 








STORM | ERROR | 




























































CAROL 16 27-0730 ~61.5 
Aug. '5h i 28-0730 -36.5 
31-1200 7 29-0730 =e5 

36 30-0730 ano 

i 31-0730 + 1,0 

EDNA 2.8 9-0730 | 10-1830 -19.5 

| Sept. '54 eee 10-0730 w1-1530 + 1.5 

11-1400 3.5 11-0730 | 11-1600 2 20 
| HAZEL 2 Ie-=1¢80.) 1520860 
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1568800 

| CONNIE 1.6 9~1530 
| Aug. 155 1.8 10-0130 
im 12=2200 hath 11-0030 
| 226 11-2100 
2.6 12-0720 | I2sB0 

DIANE 9 15-1330 | 16-1000 
Aug. '55 er 16-1330 | 16-2200 
17-1600 
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Table VIII. "DEzP WATER" FORECASTS OF TIME OF 
MAXIMUM SURGE (2 pages) 
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STORM ERROR 
time of Roo D t te (fest- 
max. Surge actual) 
GRACIE 2 26-1900 27-1900 ~/,0.0 
Sept, '59 2 27-0100 28-0100 ~34,.0 
29-1100 4 2(~0700 28-0700 ~28,0 
2 27-1300 28-1300 -22,.0 
ile 27-1900 28-11,00 -21,0 
8 28-0100 28-2000 
2 25-0700 29-0100 
1 28-1300 29-0500 
8 28-1900 29-0730 
3 29-0100 29-0830 
5 29-0700 29-1200 
0 23-0730 23-1730 
2 23-1330 23-2100 
2 23=2200 24-0200 
6 21,—-0130 24-0330 
AUDREY 3 25-0730 26-0630 
June '57 2 260-0730 26-2000 
27-1100 L, 27-0730 27-0900 













Rog - radius of 1000 mb isoline (° lat.) 

D - distance to point of max surge C lens. ) 

t — synoptic time of hurricane (DTG local) 

te - "deep water" forecast tine of max surge 
fable Vili. NDEEP WATZR" FORECASTS OF TIME OF 


MAXIMUM SURGE (Continued) 
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STORM SHALLOW | ERROR 








time of WATER t', (fcst- 
max, surge TRAVEL actuade 
a 29-0830 —51.5 
30-0930 =26.5 

31-0730 ~ 4.5 

31-2130 eo 
31-1900 + 7,0 
Za 11-0130 =e 5 
Sept. '54 11=2200 + 80 
11-1400 11-2230 + 8.5 
eae 15-0830 ae Lees 
15-1430 +2 
CONNIE 15 9-2100 = (5.0 
Aug. '55 10-1030 59.05 
12-2200 11-0630 ~29.5 
12-0230 -19.5 
12-1930 - 2,5 
DIANE LO 16-14,00 -26.0 
Aug. 55 17-0230 = 36> 

17-1600 

se -28 © 5 
-12,.0 
SOS 

Table IX. "DEEP<$SHALLOW WATER" FORECASTS 


OF TIME OF MAXIMUM SURGE (2 pages) 
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max, surge TRAVEL 
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27-2300 
2E-0500 
28-1030 
28-1700 
28-1800 
29-0000 
29-0530 
29-0900 
29-1100 
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FLOSSY 
Sept. '56 
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AUDREY 
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-20.5 
June '57 - 8,0 
27-1100 - 3 
t', - "deep-shallow water" forecast time 
of maximum surge (DTG local) 
Table IX. "DEBP~SHALLOW WATER" FORECASTS 


OF TIME OF MAXIMUM SURGE (Continued) 
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t - "surge speed = storm speed" forecast 
: time of maximum surge 
be ~ "deep water" forecast time of max. surge 
tp - "deep-shallow water" forecast time 


of maximum surge 
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values of the time that the storm would reach the coast were used, it 
is believed that the results would not be as good. The errors listed 
in Table X for the "deep water" and the "deep-shallow water" methods 
are the minimum errors which resulted, which in most cases occurred when 
the forecasts were made for hurricane positions which were quite close 
to the coast. 

The best results were obtained from letting the surge travel at 
the speed of the storm. However, it is felt desirable to have a techni-~ 
que which is independent of parameters which have to be forecast. There- 
fore a further comparison of the "deep water" and "deep-shallow water" 
methods will be discussed. 

For these techniques, Tables XI and XII show the number of cases 


falling within the error limits indicated, for forecasts made when the 


hurricanes were at different distances from the coastline. For forecasts 
made when the storm was within 300 nautical miles of the coast fairly 
reliable results were obtained. Forecasts made when the storms were 
further than 700 miles from the coast appear to have little reliability. 
Making a forecast when a hurricane is within 300 miles of the coast seems 
to leave little time for action; however, for the nine storms considered 
in this study, the least time available, from when the storm was 300 
miles distant to the occurrence of the maximum surge, was 23 hours. 
For a hurricane travelling at 30 knots, about eight hours would elapse 
between a forecast made when the storm was 300 miles at sea, and the 
time of the maximum surge. 

A comparison of Tables XI and XII show that some improvement over 
the "deep water" method is made by applying the shallow water correction. 


However, the results obtained from "deep-shallow water" forecasts are 
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MAXIMUM DISTANCE NUMBER OF CASES WITHIN 


TO STORM WHEN ABSOLUTE ERROR LIMITS INDICATED 
FORECAST IS MADE <4 hrs Lel2 12-2) ies | 2a S 
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Tab sea.0 NDBEP WATER" TIME ERRORS 
MAXIMUM DISTANCE NUMBER OF CASES WITHIN 
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FORECAST IS MADE <(penrs >2i hrs 





COO en. . 


(00, Dem. 


900 n.mi. 


Table XII. "DEXP=SHALLOW WATER" TIME ERRORS 
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not strikingly better than those derived from use of the "deep water" 
technique. Further, it is a great deal more dangerous to forecast the 
maximum surge one hour late than even several hours early, and with the 
"deep=shallow water" method many late forecasts were obtained. For 
these reasons it is not felt that the extra work involved in the "deep~ 
shallow water" technique is worth-while to the practical forecaster. 

As can be seen from the development of this forecasting technique 
there are many possible sources of error. Some sweeping assumptions 
have been made in the interest of producing a technique which requires 
a minimum of parameters, and that requires parameters which would be 
easily available to the forecaster before the arrival of the hurricane. 
However, from the results of applying the "deep water" method to the 
nine storms on which it was tested,it appears that the technique is well 
worth additional testing to determine if it can actually be used in the 


field, 
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4, SUGGESTED FUTURE RESEARCH AND CONCLUDING REMARKS 


Height Forecastin 





At present fairly reliable results can be expected from height of 
maximum surge forecasts made using the NHRP method. It appears that the 
Kajiura method would give better results than the NHRP technique if it 
were possible to determine the amplification factcr to be applied in a 
particular forecasting situation. An empirical study of this factor 
might reveal an accurate and convenient method of applying it. 

4.2 Time Forecasting. 

The method developed in this study should be tested further to 
determine its actual usefulness. With more data available, a strictly 
empirical treatment of the parameters used in this study to forecast 
the time of the maximum surge might result in a more reliable technique. 

The present measurement of the size of the storm, the radius of 
the 1000 mb isoline, can be quite inaccurate, Perhaps the substitution 
of some related parameter which maybe more correctly determined would 
improve the results. 
ho3 Data. 

In conclusion it is well to point out that the main thing that 
has held back worthwhile study in this field has been the lack of data. 
There is available a wealth of hurricane data for almost all parameters 
imaginabie except for the high waters associated with the storm. 

This problem may soon be alleviated, for by the end of 1960 the 
U. S. Weather Bureau expects to publish a series of storm surge charts. 
some of the data for this study were taken from advance copies of these 


charts, These charts show the storm track geographically and relative 


43 





to the maximum surge profiles which are drawn for a series of stations 
along the coast. Further, the Climatological Data publication has start-~ 
ed publishing similar surge data in a less elaborate form for current 
storms. 

With these additional data available we may expect great strides 
toward accurate forecasting of all aspects of the surge problem in the 


future, 


Ad 
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APPENDIX T 


THE CONNER, KRAFT, AND HARRIS HeIGHT FORECASTING ifsTHOD 


ineeprat 1957. Cee Conneryekow. deafteand D. Leemiarers [2 | 
presented a quantitative technique for forecasting the maximum height 
of the storm surge. 

In the method, which was similar to one which had been developed 
for local use by the New Orleans Weather Bureau Office several years 
previously, a single parameter was selected which would cover as many 
of the factors which affect the maximum surge as possible. A simple 
storm model was used in which it was assumed that the pressure distri- 
bution depended only on the distance from the center of the storm. In 
such a model the faximeeetind speed, V eae depends on the pressure 
deficiency as expressed by; 


Ve Ces) 


Vou 


whe e: K = constant of proportion- 
ality 
Dp. = pressure at the edge of 
n 
the storm 
Pp, = pressure ab tne center of 


the storm, 


The total setup due to the wind, h, was taken as; 
ab 
h re 


wheres B& b=constants determined 
empirically. 


These two equations were combined to obtain a suitable form for calcula= 
tion of the maximum surge height based on pressure, It had been found 
that the pressure at the edge of the storm could be considered constant 
for all storms. In this case a value of 1019 mb was selected for Ps 
and b was chosen to equal one. B was evaluated using data for Gulf of 


Mexico hurricanes and the final resulting equation took the form: 


_” joer 
| 
-—— — | ee 
‘_.——S—- --——— —b— & <——-s 
me memes ily” — oma mm mma 
eo <P aa 
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Z a0 as 
hvnax = O.154 (1019~p, } 
where: h = maximun storm surge 
max a 
height. 

Although based on Gulf coast data, this method was tested using 
Atlantic coast hurricanes and it was found satisfactory except for 
storms which travelled parallel and close to the coastline. 

A graph of the final forecasting equation was presented which 
the forecaster could enter with the lowest pressure in the hurricane 
A copy of this forecasting 


to obtain the maximum storm surge height. 


curve is included here as Figure ll. 
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Figure ll. CONNER, KRAFT, AND HARRIS 
HEIGHT FORECASTING GRAPH 


Ag 





APPENDIX II 


THE HOOVER HEIGHT FORECASTING METHOD 


Concurrently with the development of the Conner, Kraft and Harris 
study (Appendix I), Robert A. Hoover Le also developed a technique 
for forecasting the maximum height of the storm surge. 

The development of the Hoover method was essentially the same as 
for the Conner, Kraft and Harris technique, except that the analysis 
of the data was handled differently and more data were used. Hoover 
pointed out that since the actual highest tide may frequently go un- 
observed at tide stations, the regression line determined by Connor, 
Kraft, and Harris was biased toward underestimating the maximum. 
Furthermore, in the previous study many sets of hurricane data could 
not be used since the tide gages became inoperative before the maximum 
surge height was reached, 

For as many storms as data were available, Hoover constructed pro=- 
files of the surge height versus lateral distance from the storm track. 
The data was weighted in an effort to obtain profiles which would not 
reflect local effects, and where data were missing the profile was 
estimated, Upon examination of these surge profiles it was found that 
along certain parts of coastline they resembled one another more than 
along other parts of the coastline. These groups of profiles were 
e@iererore combined on a regional basis. One of these sets of surge 
profiles, that for hurricanes which entered the coasts of North Caro- 
line, South Carolina, and Georgia, is included in this study as Figure 
ee 

Hoover used these profiles to interpolate the actual highest surge 


in cases where data were missing or doubtful. From the data thus obtain- 
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Figure Ze —Surge profiles of the hurricanes which entered the coasts 
of North Carolina, South Carolina, and Georgia. a 





ed two regression lines of pressure and maximum surge height were deter- 
mined, one for the Atlantic coast and another for the Gulf coast. The 
expressions for these two regression lines were given as: 


Atlantic coast; 


hex = 0.198 \1006-p, } 
Gulf coast; 
here = 0.15! (1032~ po) 
where: h 


p 


ext = max. surge height 


lowest pressure at 


O 
center of storm 


These two regression lines, which may be used by the forecaster to 


determine the maximum surge height, are shown in the present study as 


Bieure 13, 
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Figure 13. HOOVER HEIGHT FORECASTING GRAPHS 
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APPENDIX III 


THE NATICNAL HURRICANE RESEARCH PROJECT 
INTERIM HRIGHT FORECASTING METHOD 


In August of 1959 the U. S. Weather Bureau officially published 
An Interim Hurricane Storm Surge Forecasting Guide" \é6}. The larg- 
est part of this report was devoted to the presentation of a method for 
forecasting the maximum height of the surge. 

In the development of this method it was assumed that the storm 
surge height depends on the wind speed, the fetch length, and the depth 


of water. The basic equation that was used is: 


where: h storm surge height 


va effective wind stress 
directed toward the 
shore 

L = length of the fetch 

D = depth of the water. 


It was further assumec that the effective wind stress was related 
to the pressure deficiency, (p, - Pode A constant value was assumed 
for Py the pressure at the edge of the storm. Additional data, more 
reliable than previously, were available for this study. A value of 
1025 mo. was determined for P, by a least-square fit of the data. It 
was found that making use of the length of the fetch, L, did not improve 
the results sufficiently to be considered. Several parameters were con= 
sidered for determination of the depth of the water. The most effective 
was found to be the distance from the shore to the 50~fathom curve. 


The equation which resulted from these considerations was: 


h= 0.06013 (l025-p) Se 


where: Bare central pressure of the 


storm 

S = distance between the shore 
and the 50fathom depth 
COnGour. 


a4 


This equation was plotted in the form: 


h=6 4(5,) 


where; b relative storm surge 
potential (depends on 
off-shore depths) 

a function of the cen-= 


tral pressure. 


S (p,) 


A map was provided, shown in this study as Figure 14, which showed 
values of @ for the Atlantic and Gulf coasts. A storm surge prediction 
chart was developed which gave values of surge fesene versus © and Po? 
This chart is shown here as Figure 15. 

To make a prediction of the maximum height of the storm surge the 
forecaster obtains a value of G from the Map of Relative Storm Surge 
Potential using the storm track, and then enters the Storm Surge Pre- 
diction Chart with this parameter and the lowest pressure in the storm 


to obtain the height. 


aD 





Figure 14, 
- Map showing relative storm surge potentinl 9 of various coastal 
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APPENDIX IV 


THE KAJIURA HEIGHT FORECASTING METHOD 


In December of 1959 Kinjiro Kajiura published an extensive 
study of all the phases of the storm surge problem [3]. A part of 
Kajiura's work includes a suggested technique for forecasting the maxi- 
mum height of the storm surge. 


The basic equation used in this study was: 


ae where? si = maximum storm surge height 
r = dynamic amplification 
factor 
a' = constant depending on the 
station 
P = constant density of water 
g = acceleration due to gravity 
Ap*= intensity of the storm. 


Values of a' were found that covered two rather long sections of the 
Atlantic coastline, one extending along the New England Coast southward 
to Montauk, New Yorks and the other extending from Long Island to the 
Florida peninsula. When combined with e and g these values were 0.14 
for the New England region, and 0.2 for the remainder of the Atlantic 
coast. For the Gulf coast the value of these combined constants was 
not determined but it was believed to be larger than for the Atlantic 
coast regions due to the existence of a shallow and wide shelf. 

The factor of dynamic amplification was not defined for particular 
forecasting situations by Kajiura. It was stated that r lies between 
l and 2, and for realistic situations it would seem to be around 1.5 to 
2 if the storm moved almost perpendicular to the coastline, over a wide 
continental shelf, and if the speed of storm movement was suitable. The 


results of a previous determination of the dynamic amplification factor 


58 


—_ EP 
ae ‘> > > 
- _?* -- = 
——_—_ — as == & 
a 2? eee es eS a 
a - a ny 
; ~~ ae 
- —- > 
. ee 
= ee ee + 
aa aatieiannaiio i A se 
a a | 
_ “oe 


7 7 
| anal ae 
a - a 





by Reid in 1956 for a theoretical one-dimensional model were given. 

For this one=dimensional case r was somewhere around 1,5 ~— 2.0. 

Kajiura further stated tnat the amplification factor for two dimension-~ 
al models is a little smaller than for the one-dimensional case, There 
were not enough data available at the time that this method was develop- 
ed to determine the amplification factor on an empirical basis. 


: . ae ; 
An approximate expression for Ap was taken as: 


Pate 5 (Ol ps) 


where: Bao lowest central pressure. 

This expression can be combined with the basic equation to gives 
T= ralo.s5(10l5~p,)] 
wheres: a = constant combining a’, 
e » & Zo 

From this last equation a curve (or series of curves) for the particular 
region determined by the constant, a,could be drawn if the value (or 
series of values) of r were known. Two such graphs could be used to 


forecast the maxinum height of the surge. 
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